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Abstract. Wc report on XMM-Newton observations of 
the "Black Widow pulsar", PSR B1957+20. The pulsar's 
X-ray emission is non-thermal and best modeled with a 
single powerlaw spectrum of photon index 2.031q g 6 . No 
coherent X-ray pulsations at the pulsar's spin-period could 
be detected, though a strong binary-phase dependence of 
the X-ray flux is observed for the first time. The data 
suggest that the majority of the pulsar's X-radiation is 
emitted from a small part of the binary orbit only. We 
identified this part as being near to where the radio eclipse 
takes place. This could mean that the X-rays from PSR 
B1957+20 are mostly due to intra-shock emission which is 
strongest when the pulsar wind interacts with the ablated 
material from the companion star. 



1. Introduction 

Until now, more than 1700 rotation-powered radio pulsars 
are detected. Among them are about 10% which are mil- 
lisecond pulsars (MSPs) (Manchester et al. 2005). They 
form a separate population. The majority of them re- 
sides in Globular Clusters (c.f. Bogdanov et al. 2006). 
MSPs are presumed to have been spun up in a past ac- 
cretion phase by mass and angular transfer from a bi- 
nary companion (Alpar et al. 1982). Only about one third 
of them are seen to be solitary. It is believed that they 
lost their companion, e.g. in a violent supernova event. 
All MSPs possess very short spin periods of less than 
20 ms and show a high spin stability with period deriva- 
tives in the range sa 10~ 18 — 10 -21 . MSPs are generally 
very old neutron stars with spin-down ages r = P/2P of 
~ 10 9 — 10 10 years and low surface magnetic fields in the 
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range B cx \I{PP) ~ 10 8 - 10 10 G. 

At present, about 50% of all X-ray detected rotation- 
powered pulsars are MSPs (c.f. Bogdanov et al. 2006 and 
references therein). Among them an extraordinarily rich 
astrophysics binary system which is formed by the mil- 
lisecond pulsar PSR B1957+20 and its 0.025 M Q low mass 
white dwarf companion (Fruchter, Stincbring, & Taylor 
1988). The pulsar has a spin period of 1.6 ms which is the 
third shortest among all known MSPs. Its period deriva- 



tive of P — 1.69 x 10~ 20 s s _1 implies a spin-down energy 
of E — 10 35 erg s~ l , a characteristic spin-down age of 
> 2 x 10 9 years, and a dipole surface magnetic field of 
B± = 1.4 x 10 s Gauss. The pulsar is orbiting its compan- 
ion with an orbit period of 9.16-h. Optical observations 
by Fruchter et al. (1988) and van Paradijs et al. (1988) 
revealed that the pulsar wind consisting of electromag- 
netic radiation and high-energy particles is ablating and 
evaporating its white dwarf companion star. This rarely 
observed property gave the pulsar the name black widow 
pulsar. Interestingly, the radio emission from the pulsar 
is eclipsed for approximately 10% of each orbit by mate- 
rial expelled from the white dwarf companion. For a radio 
dispersion measure inferred distance of 1.5 kpc (Taylor 
& Cordes 1993) the pulsar moves through the sky with 
a supersonic velocity of 220 km/sec. The interaction of a 
relativistic wind flowing away from the pulsar with the in- 
terstellar medium (ISM) produces an H a bow shock which 
was the first one seen around a "recycled" pulsar (Kulka- 
rni & Hester 1988). 

In 1992 Kulkarni et al. published a contour map of the 
X-ray emission of PSR B1957+20 which was derived from 
ROSAT PSPC observations. Although this ROSAT data 
were very sparse in statistics it led the authors to predict 
faint diffuse X-ray emission to be present along a cylindri- 
cal trail formed when the relativistic pulsar wind expands 
into pressure equilibrium with the interstellar medium be- 
hind the nebula. The much improved sensitivity of the 
Chandra and XMM-Newton observatories made it possi- 
ble to probe and investigate the structure and properties 
of this unique binary system in much higher detail than it 
was possible with ROSAT, ASCA or BeppoSAX. A narrow 
X-ray tail with the extent of 16 arcsec and the orientation 
to the north-east was detected from it in deep Chandra 
observations by Stappcrs et al. (2003). 

Searching in ROSAT data for a modulation of the pul- 
sar's X-ray emission as a function of its orbital phase re- 
vealed a suggestive but insignificant increases in flux be- 
fore (at phase <fi ~ 0.17) and after (at phase <f> ~ 0.4 — 0.5) 
the pulsar radio eclipse (cj> = 0.25) (Kulkarni et al. 1992). 
Taking Chandra data into account revealed a hint that 
the lowest and highest fluxes are located during and im- 
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mediately after the radio eclipse, respectively. The statis- 
tical significance of this modulation observed by Chandra, 
though, is only 98% and thus prevents any final conclusion 
on it (Stappers et al. 2003). 

2. Observations and Data Analysis 

PSR B1957+20 was observed with XMM-Newton on Oc- 
tober 31, 2004 for a 30 ksec effective exposure. In this ob- 
servation, the EPIC-MOS1 and MOS2 instruments were 
operated in full-frame mode using the thin filter to block 
optical stray light. The EPIC-PN detector was setup to 
operate in the fast timing mode. Because of the reduced 
spatial information provided by the PN in timing mode we 
use the MOS1/2 data for imaging and spectral analysis of 
the pulsar and its diffuse X-ray nebula while the PN data 
having a temporal resolution of 0.02956 ms allowed us to 
search for X-ray pulsations from the pulsar. All the data 
were processed with the XMM-Newton Science Analysis 
Software (SAS) package (Version 6.5.0). Spatial and spec- 
tral analyses were restricted to the 0.3 — 10.0 kcV energy 
band while for timing analysis events were selected for the 
energy range 0.3 — 3.0 keV. 

2.1. Spatial Analysis 

Figure 1 shows the combined EPIC-MOS1/2 image of the 
PSR B1957+20 system. The image was created with a bin- 
ning factor of 6 arcsec and by using an adaptive smoothing 
algorithm with a Gaussian kernel of a < 4 pixels in order 
to better make visible faint diffuse emission. The extent of 
the diffuse emission with its orientation to the north-east 
is about 16 arcsec which is consistent with the previous 
result derived from the Chandra Observations. However, 
the detailed structure of the X-ray emission from XMM- 
Newton can not be as clearly seen as from Chandra due 
to the 10 times wider Point Spread Function (PSF) of 
XMM-Newton. 

Inspecting the XMM-Newton MOS1/2 image two faint 
features (denoted as A and B) which contribute to only 
about 3% of the total X-ray flux are apparent. In order 
to investigate whether this faint features are associated 
with nearby stars we inspected the Digitized Sky Survey 
data (DSS) and the USNO-B1.0 Catalogue for possible 
sources. These catalogues which are limited down to 22 
mag (Krongold, Dultzin-Hacyan & Marziani 2001) and 21 
mag (Monet et al. 2003) do not reveal possible counter- 
parts. These features are not seen in the Chandra image 
though (Stappers et al. 2003). 

2.2. Spectral Analysis 

Combined EPIC-MOS1/2 data of PSR B1957+20 were 
extracted from a circle of 30 arcsec radius centered at 
the pulsar position (RA (J2000) = 19 /l 59" l 36 s .77, Dec 
= 20°48'15" .12). The selection region contains about 85% 




Fig. 1. XMM-Newton MOS1/2 image of the 
PSR B1957+20 -WD system with contour lines 
overlaid. The contour lines arc at the levels of 
(5.1,6.1,9.4,16.3,32.4,69.0) x 10~ 6 cts s-^rcscc" 2 . 
The position of the pulsar is indicated. 

of all source counts. Background photons were selected 
from a source-free region near to the pulsar position. Re- 
sponse files were derived by using the XMM-Newton SAS 
tasks RMFGEN and ARFGEN. 

After subtracting background photons, in total 338 
sources counts were available for a spectral analysis. The 
extracted spectra were binned with at least 30 source 
counts per bin. Assuming that the emission originates 
from the interaction of the pulsar wind with the ISM or 
with the stellar wind we expect synchrotron radiation to 
be the emission mechanism of the detected X-rays. To test 
this hypothesis we fitted the spectrum with a power law 
model. Indeed, this model describes the observed spec- 
trum with a reduced \l of 1-09 (for 8 D.O.F.). The 
photon-index is found to be a = 2. 03^36 ■ The column 
absorption Nh is 8.0 x 10 20 cm~ 2 . For the normalization 
at 1 keV we find 1.5t°J x 10 -5 photons keV~ 1 cm~ 2 sec _1 
(l-<7 confidence for 1 parameter of interest). The spectrum 
and the fit residuals are shown in Figure 2. 

The unabsorbed X-ray fluxed derived from the best fit- 
ting model parameters is f x = 8.35 x 10~ 14 erg s _1 cmT 2 
and f x = 7.87 x 10~ 14 erg s" 1 cjti- 2 in the 0.3 - 10 
keV and 0.1 — 2.4 keV energy band, respectively. The X- 
ray luminosities in these energy bands - calculated for 
a pulsar distance of 1.5 kpc - are L x (0.3 — lO.OkeV) = 
2.24 x 10 31 erg s" 1 and L x (0.l - 2.4kcV) = 2.12 x 
10 31 erg s , respectively. The conversion efficiency L x /E 
in the 0.1 - 2.4 kcV band is found to be ~ 2.12 x 10~ 4 . 

2.3. Timing Analysis 

The EPIC-PN camera observed the pulsar in the fast tim- 
ing mode. In this mode the spatial and spectral informa- 
tion from a 64 x 199 CCD pixel array is condensed into 
a one dimensional 64 x 1 pixel array (ID-image), i.e. the 
spatial information in Y-direction is lost due to the con- 
tinuous read-out of the CCD. The complete photon flux 
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channel energy (keV) 

Fig. 2. Energy spectrum of PSR B1957+20 obtained from 
the XMM-Newton MOS1/2 data. The plot shows the X- 
ray spectrum fitted with an absorbed power law model 
(upper panel) and contribution to the \ 2 fit statistic (lower 
panel) . 



(source plus DC emission from foreground or background 
sources located along the read-out direction) is accumu- 
lated and collapsed in the final ID-image, severely reduc- 
ing the signal-to-noise ratio of pulsed emission and pre- 
venting the detection of weak X-ray pulsations from the 
target of interest. 

In order to search for X-ray pulsations from PSR 
B1957+20 we extracted 1639 counts from the CCD 
columns 33 — 41 in which the pulsar got located. In or- 
der to increase the signal-to-noise ratio we restricted the 
analysis to the energy range 0.3 — 3.0 keV. Below and be- 
yond this energy band the accumulative sky and instru- 
ment background noise exceeds the contribution from the 
pulsar itself (c.f. Becker & Aschenbach 2002). However, 
still 80 % of the counts are estimated to be derived from 
the background. 

The photon arrival times were corrected to the so- 
lar system barycentre with the BARYCEN tool (version: 
1.17.3, JPL DE200) of the SAS package. As the pulsar is in 
a binary we corrected for the orbital motion of the pulsar 
by using the method of Blandford & Teukolsky (1976). 

As millisecond pulsars are known to be extremely sta- 
ble clocks we used the pulsar ephemeris from the ATNF 
Catalogue, / = 622.122030511927 Hz and / = -6.5221 x 
10~ 15 sec- 2 (at MJD = 48196.0) to perform a period fold- 
ing. By using the Z\ statistics (Buccheri et al. 1983) with 
the harmonics number (n) from one to ten no significant 
signal was detected at the radio spin period extrapolated 
for the epoch of the XMM-Newton observation. Restrict- 
ing the period search to the various smaller energy bands 
did not change the result. A pulsed fraction upper limit 
of 9% (1-er) is deduced by assuming a sinusoidal pulsed 
profile. 

Arons & Tavani (1993) predicted that depending on 
the flow speed and the degree of absorption and/or scat- 
tering by the companion wind the X-ray emission from 
PSR B1957+20 increases by up to a factor of 2.2 at the 




Orbital Phase 



Fig. 3. X-ray emission from PSR B1957+20 within 0.3 - 
3.0 keV as function of the pulsar's orbital phase (<fi). = 
corresponds to the ascending node of the pulsar orbit. The 
upper curve was obtained from the XMM-Newton EPIC- 
PN (background level at 88 cts/bin). The lower lightcurve 
is obtained from the MOS1/2 data. Phase bins with zero 
counts correspond to phase angles not covered in the ob- 
servation. 



orbital phase before and after the radio eclipse. In order to 
test this prediction we created a lightcurve by binning all 
events in bins of 2.2 ksec width. With an effective exposure 
time of about 30 ksec and an orbit period of 9.16-h the 
XMM-Newton data cover roughly one binary orbit. The 
lightcurves resulting from the EPIC-PN and MOS1/2 data 
are shown in Figure 3. For the lightcurve deduced from the 
PN-data (upper solid curve) it is clearly seen that before 
the radio eclipse, i.e. between the orbital phase 0.05 and 
0.25, the X-ray emission increases. The emission in the 
highest bin is about a factor of 2.5 stronger than at other 
orbital phase angles before and after the radio eclipse. 
This is in agreement with the predictions made by Arons 
& Tavani (1993). 

Although it is the first time that a significant X-ray 
flux modulation from PSR B1957+20 is observed, the flux 
increase near the phase of the radio eclipse is virtually 
only seen in the PN data. Owing to the short observation 
time of 30 ksec, which is less than the time of one full 
binary orbit, the orbit angles 0.30 — 0.38 and 0.18 — 0.36 
are not covered at all. It thus turns out that it was a 
great fortune that the EPIC-PN, more or less by chance, 
covered the orbital phase range of the radio eclipse and 
provided us evidence for the flux enhancement while the 
MOS1/2 detectors were already switched off. Table 1 lists 
the first and the last photon arrival times recorded by the 
MOS1/2 and EPIC-PN detectors and the corresponding 
pulsar binary orbital phase. 



3. Summary and Discussion 

Interaction between relativistic pulsar winds which carry 
away the rotational energy of pulsars and the surrounding 
medium is expected to create detectable X-ray emission. 
Indeed, there are about 30 pulsar wind nebulae (PWNe) 
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Table 1. List of the first and the last photon arrival times and its corresponding orbital phase of PSR B1957+20 
for the MOS1/2 and PN detectors. Arrival times are corrected to the solar system barycentre. The orbital phase is 
measured from the time of ascending node. 



Date set 


1st photon 


last photon 


Duration (s) 




Time (MJD) Phase 


Time (MJD) Phase 




MOS1 


53309.9688 0.3605 


53310.2830 0.1830 


27143.4 


MOS2 


53309.9666 0.3548 


53310.2846 0.1871 


27467.1 


PN 


53309.9791 0.3873 


53310.3301 0.3063 


30327.4 



currently detected in the X-ray band (e.g. Kaspi, Roberts 
& Harding 2006, Gaensler & Slane 2006, Kargaltsev & 
Pavlov 2006). However, these PWNc arc all powered by 
young and powerful pulsars with spin-down energies of 
more than ~ 3.6 x 10 36 erg s _1 . Until now, only two 
MSPs are known to have X-ray nebulae. They are PSR 
B1957+20 (Stappers et al. 2003) and PSR J2124-3358 
(Hui & Becker 2006). Both of them have tail-like struc- 
tures behind the moving pulsars. These trails could be 
associated with the shocked relativistic wind confined by 
the ram pressure of the ambient ISM. 

The XMM-Newton data of PSR B1957+20 have pro- 
vided observational evidence for a strong dependence of 
the pulsar's X-ray emission on its binary orbital phase. It 
is the first time that a significant X-ray flux modulation 
from PSR B1957+20 is observed. The emission near to the 
radio eclipse is supposed to be beamed in a forward cone 
because the shocked fluid is accelerated by the pressure 
gradient as it flows around the eclipse region. Relativistic 
beaming would tend to give the maximum flux just be- 
fore and after the pulsar eclipse which for PSR B1957+20 
is at orbital phase = 0.25. Arons & Tavani (1993) pre- 
dicted that the immediately downstream flow velocity of 
the shocked pulsar wind from the shock area along the line 
between the pulsar and the ablating companion is about 
c/3 and is even higher behind the relativistic shock when 
it passes around the companion. The tenuous relativis- 
tic plasma accelerates as it flows around the companion, 
possibly passing through a sonic transition to leave the 
binary system with the velocity larger than c/v3. Due 
to Doppler boosting, the probable post-shock velocities of 
the relativistic wind then suggest X-ray emission variation 
around the orbital phase by a factor between 1.3 (v = c/3) 
and 2.2 (v = c/\fi). This numbers were estimated with- 
out considering the effect of absorption and/or scattering 
within the binary. On the contrary, the X-ray emission at 
the eclipse may be reduced because of the obscuration of 
the shock by the companion. 

However, given the limited photon statistics of the 
XMM-Newton data it is not possible to investigate any 
binary-phase resolved imaging and spectral variation as 
a function of orbit phase or to determine the exact ge- 
ometry of the peak emission. Such analysis would allow 
us to investigate whether the X-ray emission from PSR 
B1957+20 is present during all orbit angles or virtually 



only near to the radio eclipse while diffuse X-ray emission 
from the PWN is present at any orbital angle. Indeed, the 
later is suggested by the current data and a confirmation 
would have a severe impact on our understanding of the 
pulsar's X-ray emission properties. As the present XMM 
observation covers barely one binary orbit, we stress that 
it can not be fully excluded that the increase in photon 
flux near the orbital angle of the radio eclipse is due to 
a single "burst like" event or that the peak flux emis- 
sion varies from orbit to orbit or on longer time scales. A 
repeated coverage of the binary orbit in a longer XMM- 
Newton observation and a comparison with the 2004 data 
would answer this question immediately. In addition, this 
would provide us not only a better photon statistic but 
also would allow us to determine the emission geometry 
with a much higher accuracy than currently possible. 
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